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We probe the local detection efficiency in a nanowire superconducting single-photon detector along 
the cross-section of the wire with a spatial resolution of 10 nm. We experimentally find a strong 
variation in the local detection efficiency of the device. We demonstrate that this effect explains 
previously observed variations in NbN detector efficiency as function of device geometry. 


Nanowire superconducting single-photon detectors 
(SSPDs) consist of a superconducting wire of nanoscale 
cross-section [1], typically 4 nm by 100 nm. Photon de¬ 
tection occurs when a single quantum of light is absorbed 
and triggers a transition from the superconducting to the 
normal state. SSPDs have high efficiency, low jitter, low 
dark count rate and fast reset time [2], and are therefore 
a key technology for, among others, quantum key distri¬ 
bution [3], interplanetary communication [4] and cancer 
research [5]. 

Although progress has been made recently, the under¬ 
lying physical mechanism responsible for photon detec¬ 
tion on the nanoscale is still under active investigation. 
A combination of theory El El, experiments EHini and 
simulations [TT1[T2] on NbN SSPDs indicates that the ab¬ 
sorption of a photon destroys Cooper pairs in the super¬ 
conductor and creates a localized cloud of quasiparticles 
that diverts current across the wire. This makes the wire 
susceptible to the entry of a superconducting vortex from 
the edge of the wire. Energy dissipation by this moving 
vortex drives the system to the normal state. 

An important and unexpected implication of this de¬ 
tection model is a nanoscale position variation in the 
photodetection properties of the device. The conditions 
at the entry point of the vortex determine the energy 
required for it to cross the wire. This causes photons 
absorbed close to the edge to have a local detection effi¬ 
ciency (LDE) [13] compared to photons absorbed in the 
center of the wire [T2|. This effect has practical impli¬ 
cations for the operation of SSPDs, since it represents 
a potential limitation of the detection efficiency. In ad¬ 
dition, SSPDs have been proposed for nanoscale sensing, 
either in a near-held optical microscope conhguration El 
or as a subwavelength multiphoton probe El, where this 
effect would be of major importance for the properties of 
such a microscope. While this effect has been predicted 
theoretically, clear experimental evidence is missing. 
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In this work, we experimentally explore the nanoscale 
variations in the intrinsic response of the detector. We 
spatially resolve the LDE with a resolution of approxi¬ 
mately 10 nm, better than A/50, using far-held illumi¬ 
nation only. We hnd that our results are qualitatively 
consistent with numerical simulations E1II2- Our re¬ 
sults provide excellent quantitative agreement with ex¬ 
periments that indicate a polarization dependence in the 
LDE that was hitherto not understood El- 

The key technique used in this work is a differential 
polarization measurement that probes the IDE of the de¬ 
tector (see Eigure I). The technique is based on the fact 
that polarized light is absorbed preferentially in differ¬ 
ent positions for the two orthogonal polarizations, due to 
differences in boundary conditions. Using this technique, 
we achieve selective illumination of either the sides or the 
middle of the wire, which we use to probe the intrinsic 
photodetection properties of our device on the nanoscale. 

However, it is well known that changing the polar¬ 
ization results in a change in overall optical absorption 
probability [T6H2Q] . In order to correct for this effect, we 
must separate the probability that a photon is absorbed 
from the probability that an absorbed photon causes a 
detection event. To separate changes in optical absorp¬ 
tion from the intrinsic LDE effects which are presently 
of interest, we use quantum detector tomography (QDT) 
[Hl[i|2lH29]. 

QDT records the detector response to a set of known 
quantum states of light and distills from these measure¬ 
ments the detection probability for different photon num¬ 
bers. As we showed previously [28|, this procedure al¬ 
lows us to unambiguously separate the intrinisic one- 
photon detection probability pi from the probability r] 
that a photon is absorbed m- We found that r] is al¬ 
most independent of detector bias current and that the 
value is consistent with the geometric area of the de¬ 
tector. Hence, we identify pi with the detection prob¬ 
ability conditional on photon absorption, which we dub 
the internal detection efficiency (IDE). By construction, 
IDE = f LDE{x)A{x)dx/ f A(x)dx, where the integral 
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runs along the cross-section of the wire and A{x) repre¬ 
sents the absorption probability density. 
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Figure 1. Sketch of our experiment. Top: Absorption as a 
function of position in the wire for parallel (||) and perpendic¬ 
ular (J_) polarizations, calculated with an FDTD method (see 
text). Inset: Sketch showing the two polarizations. The red 
arrow represents the polarization of the electric field. Bottom: 
Experimental setup. Our laser pulses are tuned in intensity 
by a variable attenuator consisting of two crossed polarizers 
and a A/2 wave plate. Polarization is set by an additional 
A/2 wave plate. The image is a SEM micrograph of a detec¬ 
tor nominally identical to the one used in this experiment. 
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Eigure 2. Experimental results on the polarization depen¬ 
dence of the internal detection efficiency (IDE). Each graph 
represents one data point in the graph below it. Top panel, 
blaek squares: IDE as a function of polarization at an excita¬ 
tion wavelength of 1500 nm at h = 25.8 jaA. The red line is 
a sine fit. Top panel, blue eireles: Absorption probability as 
a function of polarization. Middle panel: IDE as a function 
of bias current. We find that the IDE is independent of bias 
current. The black line represents the weighted average of 
our measurements at different bias currents. Bottom panel: 
IDE as a function of illumination wavelength. We show the fit 
of the observed polarization-dependent IDE to the position- 
dependent local detection efficiency (LDE). 


We perform our experiments on a 100 nm long, 150 
nm wide NbN bridge patterned from a 5 nm-thick NbN 
film sputtered on a GaAs substrate [31] (G = 28 /iA). 
We apply a bias current and read out the detector with a 
bias-tee to separate high-frequency detection pulses from 
the DC bias current. The resulting pulses are fed to a 
series of RF amplifiers and a pulse counter. At each com¬ 
bination of bias current, photon energy and polarization, 
we record the detector count rate as a function of input 
intensity. The probe states were prepared by a broad¬ 
band pulsed laser (Fianium, repetition rate 20 MHz) out 
of which we select a narrow wavelength band with di¬ 
electric filters [9|. We prepare the desired intensity and 
polarization by first attenuating the light with a com¬ 
bination of two crossed polarizers and a half-wave plate, 
and then setting the polarization with an additional wave 
plate (see Figure 1) [3^ . 

Figure 2 shows our measured IDE and 77 , as function 
of polarization and wavelength. The top panel shows 
the IDE of our device as a function of the polarization 
of light with A = 1500 nm wavelength. Our experiments 
show that the IDE and absorption probability oscillate in 
phase when the polarization is rotated, with a minimum 
at perpendicular polarization and a maximum at parallel 
polarization. This demonstrates that absorption of TM- 
polarized photons is less likely to result in a detection. 
This polarization is absorbed preferentially in the middle 


of the wire. Our result therefore confirms earlier, prelim¬ 
inary results HSIIMIIM] that suggested that the edges of 
the detector are more efficiently photodetecting than the 
center of the wire. The error bars in this panel correspond 
to the standard deviation of a series of independent ex¬ 
periments. The middle panel of Figure 2 shows that the 
visibility of the IDE Vide = {Pmax-Pmin)l{Pmax +Pmin) 
is independent of bias current. We can therefore as¬ 
sociate one visibility to a particular illumination wave¬ 
length. For A = 1500 nm we find Vide = 0.09. The 
bottom panel shows the wavelength dependence of the 
visibility. We find that longer wavelengths have higher 
visibility in the IDE. In this panel, we show a fit to our 
position-dependent LDE, which we will discuss below. 

In the second part of this work, we will discuss our 
reconstruction of the LDE profile from these measure¬ 
ments. We make use of the fact that the IDE is given by 
the LDE multiplied by the optical absorption probability, 
integrated across the wire. Our strategy is to take the 
absorption profiles as given - since they are well studied 
- and to take the LDE profile as a free parameter and fit 
it to our experimental data. 

To calculate the optical absorption distributions, we 
perform a series of numerical simulations at differ¬ 
ent wavelengths using a finite-difference time domain 
(FDTD) method (RSOFT Fullwave). We perform a 2D 
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Figure 3. Local detection efficiency (LDE) and threshold cur¬ 
rent in an SSPD, for A = 1500 nm. Top panel: observed 
positional variations in the LDE of an SSPD, as a function 
of bias current. Inset: IDE, as a function of bias current. 
Bottom panel: Threshold current Ith- The red curve shows 
the experimentally observed value, the black curve shows the 
theoretical value obtained from numerical modelling of the 
detection event. 


simulation, and compute the electric field in a 150 nm 
wide, 5 nm thick NbN wire on a semi-infinite GaAs sub¬ 
strate and an 80 nm thick HSQ layer on top of the NbN 
wire. The refractive index of NbN deposited film on 
GaAs is derived from spectroscopic ellipsometry mea¬ 
surements [35] jSnj. In Figure 1, the result of this cal¬ 
culation is shown for A = 1500 nm. 

In order to combine information from different wave¬ 
lengths - which is required to find the LDE - we must 
posit some relation between photon energy and detec¬ 
tion probability m Our ansatz is motivated by the 
experimental observation that for low detection proba¬ 
bilities, the detection probability depends exponentially 
on bias current, by our earlier work on the energy-current 
relation in SSPDs, and on the numerical simulations de¬ 
scribed below EEHnilll]. We model the bias current 
dependence of the LDE as LDE(x) = min{l,exp (/5 — 
Ith{x))/I*}, with Ith{x) = Ic- l'{x)E, where h is the 
applied bias current, G is the critical current, E is the 
photon energy, and E = 0.65 jaA is an experimentally de¬ 
termined current scale, ^'{x) is the local energy-current 
interchange ratio, which parametrizes the detection prob¬ 
ability of the wire at different excitation wavelengths. 
These simple, empirical expressions enable us to com¬ 
bine our wavelength- and, polarization-dependent IDE 
measurements, and compute the LDE. 

The top panel of Figure 3 shows the resulting LDE 
profile for different bias currents, for a wavelength of 
1500 nm. We obtained this result by fitting ^{x) to the 


experimentally determined visibility of the polarization- 
dependent IDE in the range A = 500 — 1500 nm, using 
the calculated absorption profiles. We find that the LDE 
has a high value at the edges of the wire, up to a point 
roughly 40 nm from the edge. From there, the detection 
efficiency decreases; it is reduced by two orders of mag¬ 
nitude at the center of the wire. The LDE is current- 
dependent, but saturation sets in around A = 25.5 jaA 
{h/Ic = 0.91). 

The inset in the top panel of Figure 3 shows the effect 
of this saturation on the IDE. There are three regimes: 
a rolloff regime, where the detection property depends 
exponentially on bias current, a plateau regime at high 
currents, where LDE = 1, implying IDE = 1, and an 
intermediate regime of slowly increasing detection prob¬ 
ability. These regimes are marked by dashed lines in 
the inset. In the middle regime, parts of the detector 
are fully photodetecting, while other parts are still in a 
fluctuation-assisted regime (38). The variations between 
the experimental data and the values calculated from our 
fit are less than a factor of 2, which demonstrates the 
self-consistency of our results. We have also confirmed 
that the observed LDE profile reproduces the fact (see 
Figure 2, middle panel) that in our measurement range, 
the visibility of the polarization-dependent IDE does not 
depend on the applied bias current. From these observa¬ 
tions, we conclude that our description is able to explain 
all of the available experimental data regarding polariza¬ 
tion and bias current dependence of the photoresponse 
of this detector. 

We estimate the resolution in our measurement by 
varying the distance Ax at which we specify j'{x). In 
particular, we find that for Ax ~ rc/2, we are unable to 
explain our experimentally observed data. Therefore, we 
conclude that the LDE profile varies with a length scale 
smaller than the dimensions of the wire. We find that 
we achieve the best fit at Ax = 10 nm. From this, we 
conclude that this is the resolution of our experiment. 

The bottom panel of Figure 3 shows the threshold cur¬ 
rent Ith{x). The red curve shows the experimental value, 
from which the LDE(x) shown in the top panel is de¬ 
rived. The black curve in the bottom panel of Figure 3 
shows an independent, ab initio calculation of the thresh¬ 
old current based on the model described below. This 
ab initio calculation of the position-dependent detection 
probability is based on a numerical model mm that 
determines the threshold current for the detection of an 
absorbed photon of a given wavelength using a combina¬ 
tion of quasiparticle diffusion, current displacement and 
vortex entry [39] . We find reasonable agreement between 
our observed experimental results and the theoretical val¬ 
ues. From this, we conclude that this model captures the 
essential physics of the detection process in SSPDs. 

Our ab initio calculation gives a physical explanation 
for the enhanced efficiency at the edges of the wire in 
terms of our microscopic model. Comparing a photon 
absorption in the center of the wire to one at the edge, 
there are two differences. First, for an absorption event 
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Figure 4. Internal detection efficiency (IDE) of a set of SSPD 
meander devices. The solid symbols represent our calculation, 
and the open symbols represent experimental data from HD. 
Red symbols represent perpendicular (T) polarization, and 
black symbols represent parallel (||) polarization. The shape 
of the symbol represents the wire and pitch of the detector (see 
legend). The error bars on the experimental data represent 
the spread in properties between detectors of the same design. 
The diagonal line represents the case IDE = 1. Inset Ratio 
of overall absorption for the two polarizations. 


at the edge, the current density at the edge of the wire 
is reduced, due to the reduction in the number of su¬ 
perconducting electrons However, this is more than 
compensated by the reduction of the vortex self-energy, 
which is proportional to Vortices enter more easily 
when the superconductivity is weakened at their entry 
point, and that makes the detector more efficient at the 
edges. 

We note that there is some disagreement in theoret¬ 
ical literature about the predicted shape of the LDE 
curve. The alternative model of Zotova et al. m, 
which is based on the Ginzburg-Landau formalism, nat¬ 
urally takes into account vortex entry. However, it dis¬ 
regards quasiparticle diffusion and implements a hotspot 
with hard boundaries. The results from this model dis¬ 
agree strongly with our experimental results: there, a W- 
shaped threshold current profile is predicted, with thresh¬ 
old currents at the edges almost as high as in the center 
of the wire. The discrepancy between the models occurs 
precisely at the point where their ’hard’ hotspot touches 
the edge of the wire. We speculate that both models, if 
refined more, will likely converge. 

In the final part of this work, we discuss the implica¬ 
tion of these results to the construction and operation 
of SSPD devices. In the work of Anant et al. m, the 


IDE of a series of meander SSPDs of different wire width 
and pitch (wire separation) was measured. Using similar 
methods as described above, we compute the position- 
dependent optical absorption of such structures for the 
film properties given in that work. We use the observed 
LDE curve from our experiment, and compute the IDE 
and overall detection probability of such devices, using 
the same expressions which we used for our sample to 
obtain the LDE from the IDE and the optical absorption 

El- 

Figure 4 shows that this prediction agrees with the 
data. The calculation confirms the general claim made 
by Anant et al. that parallel polarization has a higher 
IDE than perpendicular. Beyond that, we are also able 
to compute the IDE for each device independently. The 
inset of Figure 4 shows the ratio of the overall efficiency 
for the two polarizations, which factors out the optical 
absorption as well as - to first order - the effect of bias 
current. It is therefore the most direct test of our IDE. 
We find excellent agreement. This demonstrates that we 
have achieved quantitative understanding of the inter¬ 
nal properties of an SSPD. It also demonstrates that our 
results are neither limited to the SSPD on which we per¬ 
formed our experiment nor by our single-wire geometry. 

To obtain a nonunity IDE, we must assume that the 
highly efficient detectors reported in m were not biased 
to their critical current. It is well known that the pres¬ 
ence of current crowding in the bends of a wire can cause 
reduction of the device critical current by as much as 40% 
[42], with typical values of 10-20% [43]. We have assumed 
Ib/Ic ~ 0-9 for all devices to produce Figure 4. This 
demonstrates that our results are relevant for the kind 
of SSPDs which are used for applications, at the typical 
currents at which they are operated. It also shows that 
our results can be used as an all-optical method for mea¬ 
suring the amount of current crowding in these devices. 
Since the differences between the two optical absorption 
profiles become larger at longer wavelengths, we expect 
our results to be particularly relevant for the engineering 
of SSPDs at mid-infrared wavelengths. 

In conclusion, we have demonstrated that the local de¬ 
tection efficiency of an SSPD depends on the position 
along the cross-section at which the photon is absorbed. 
We have probed this effect with a resolution of approxi¬ 
mately 10 nm, and found agreement with theoretical cal¬ 
culations done in the context of the diffusion-based vortex 
crossing model. From this, we conclude that this model 
contains the essential features for a complete microscopic 
picture of the detection model in SSPDs. We have com¬ 
pared these predictions of our work to results reported on 
devices used for applications and found good agreement, 
demonstrating the relevance of our results for SSPD en¬ 
gineering. These results enable quantitative modeling of 
the internal properties of SSPDs. 
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